Experimental approaches to assess melanocytes mosaicism in segmental vitiligo

Abstract Vitiligo is an autoimmune disease of the skin that results in localized or disseminated white macules. One common feature of several existing classification protocols is the distribution of the disease into two main subtypes, non-segmental vitiligo (NSV) and segmental vitiligo (SV). SV is characterized by depigmentation spreading within one or more skin segments while NSV is widespread. Several clinical-epidemiological observations suggest that SV has distinct autoimmune pathophysiology compared to NSV. Furthermore, the clinical distribution pattern of SV lesions closely resembles other melanocyte mosaicism diseases. These observations led us to hypothesize that SV is caused by a localized autoimmune reaction targeting epidermal mosaicism melanocytes. Here, we proposed examples of experimental approaches to assess mosaicism in SV patients.

Allele-dependent interaction of LRRK2 and NOD2 in leprosy.

Leprosy, caused by Mycobacterium leprae, rarely affects children younger than 5 years. Here, we studied a multiplex leprosy family that included monozygotic twins aged 22 months suffering from paucibacillary leprosy. Whole genome sequencing identified three amino acid mutations previously associated with Crohn's disease and Parkinson's disease as candidate variants for early onset leprosy: LRRK2 N551K, R1398H and NOD2 R702W. In genome-edited macrophages, we demonstrated that cells expressing the LRRK2 mutations displayed reduced apoptosis activity following mycobacterial challenge independently of NOD2. However, employing co-immunoprecipitation and confocal microscopy we showed that LRRK2 and NOD2 proteins interacted in RAW cells and monocyte-derived macrophages, and that this interaction was substantially reduced for the NOD2 R702W mutation. Moreover, we observed a joint effect of LRRK2 and NOD2 variants on Bacillus Calmette-Guérin (BCG)-induced respiratory burst, NF-κB activation and cytokine/chemokine secretion with a strong impact for the genotypes found in the twins consistent with a role of the identified mutations in the development of early onset leprosy.

Experimental approaches to assess melanocytes mosaicism in segmental vitiligo.

Vitiligo is an autoimmune disease of the skin that results in localized or disseminated white macules. One common feature of several existing classification protocols is the distribution of the disease into two main subtypes, non-segmental vitiligo (NSV) and segmental vitiligo (SV). SV is characterized by depigmentation spreading within one or more skin segments while NSV is widespread. Several clinical-epidemiological observations suggest that SV has distinct autoimmune pathophysiology compared to NSV. Furthermore, the clinical distribution pattern of SV lesions closely resembles other melanocyte mosaicism diseases. These observations led us to hypothesize that SV is caused by a localized autoimmune reaction targeting epidermal mosaicism melanocytes. Here, we proposed examples of experimental approaches to assess mosaicism in SV patients.

A systems biology approach identifies candidate drugs to reduce mortality in severely ill patients with COVID-19.

Despite the availability of highly efficacious vaccines, coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lacks effective drug treatment, which results in a high rate of mortality. To address this therapeutic shortcoming, we applied a systems biology approach to the study of patients hospitalized with severe COVID. We show that, at the time of hospital admission, patients who were equivalent on the clinical ordinal scale displayed significant differential monocyte epigenetic and transcriptomic attributes between those who would survive and those who would succumb to COVID-19. We identified messenger RNA metabolism, RNA splicing, and interferon signaling pathways as key host responses overactivated by patients who would not survive. Those pathways are prime drug targets to reduce mortality of critically ill patients with COVID-19, leading us to identify tacrolimus, zotatifin, and nintedanib as three strong candidates for treatment of severely ill patients at the time of hospital admission.

Deep resequencing identifies candidate functional genes in leprosy GWAS loci.

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.

Alveolar macrophages from persons living with HIV show impaired epigenetic response to Mycobacterium tuberculosis.

Persons living with HIV (PLWH) are at increased risk of tuberculosis (TB). HIV-associated TB is often the result of recent infection with Mycobacterium tuberculosis (M. tuberculosis) followed by rapid progression to disease. Alveolar macrophages (AMs) are the first cells of the innate immune system that engage M. tuberculosis, but how HIV and antiretroviral therapy (ART) affect the anti-mycobacterial response of AMs is not known. To investigate the impact of HIV and ART on the transcriptomic and epigenetic response of AMs to M. tuberculosis, we obtained AMs by bronchoalveolar lavage from 20 PLWH receiving ART, 16 control subjects who were HIV-free (HC), and 14 subjects who received ART as preexposure prophylaxis (PrEP) to prevent HIV infection. Following in vitro challenge with M. tuberculosis, AMs from each group displayed overlapping but distinct profiles of significantly up- and downregulated genes in response to M. tuberculosis. Comparatively, AMs isolated from both PLWH and PrEP subjects presented a substantially weaker transcriptional response. In addition, AMs from HC subjects challenged with M. tuberculosis responded with pronounced chromatin accessibility changes while AMs obtained from PLWH and PrEP subjects displayed no significant changes in their chromatin state. Collectively, these results revealed a stronger adverse effect of ART than HIV on the epigenetic landscape and transcriptional responsiveness of AMs.

Underwhelming or Misunderstood? Genetic Variability of Pattern Recognition Receptors in Immune Responses and Resistance to Mycobacterium tuberculosis.

Human genetic control is thought to affect a considerable part of the outcome of infection with Mycobacterium tuberculosis (Mtb). Most of us deal with the pathogen by containment (associated with clinical "latency") or sterilization, but tragically millions each year do not. After decades of studies on host genetic susceptibility to Mtb infection, genetic variation has been discovered to play a role in tuberculous immunoreactivity and tuberculosis (TB) disease. Genes encoding pattern recognition receptors (PRRs) enable a consistent, molecularly direct interaction between humans and Mtb which suggests the potential for co-evolution. In this review, we explore the roles ascribed to PRRs during Mtb infection and ask whether such a longstanding and intimate interface between our immune system and this pathogen plays a critical role in determining the outcome of Mtb infection. The scientific evidence to date suggests that PRR variation is clearly implicated in altered immunity to Mtb but has a more subtle role in limiting the pathogen and pathogenesis. In contrast to 'effectors' like IFN-γ, IL-12, Nitric Oxide and TNF that are critical for Mtb control, 'sensors' like PRRs are less critical for the outcome of Mtb infection. This is potentially due to redundancy of the numerous PRRs in the innate arsenal, such that Mtb rarely goes unnoticed. Genetic association studies investigating PRRs during Mtb infection should therefore be designed to investigate endophenotypes of infection - such as immunological or clinical variation - rather than just TB disease, if we hope to understand the molecular interface between innate immunity and Mtb.

The complex pattern of genetic associations of leprosy with HLA class I and class II alleles can be reduced to four amino acid positions.

Leprosy is a chronic disease caused by Mycobacterium leprae. Worldwide, more than 200,000 new patients are affected by leprosy annually, making it the second most common mycobacterial disease after tuberculosis. The MHC/HLA region has been consistently identified as carrying major leprosy susceptibility variants in different populations at times with inconsistent results. To establish the unambiguous molecular identity of classical HLA class I and class II leprosy susceptibility factors, we applied next-generation sequencing to genotype with high-resolution 11 HLA class I and class II genes in 1,155 individuals from a Vietnamese leprosy case-control sample. HLA alleles belonging to an extended haplotype from HLA-A to HLA-DPB1 were associated with risk to leprosy. This susceptibility signal could be reduced to the HLA-DRB1*10:01~ HLA-DQA1*01:05 alleles which were in complete linkage disequilibrium (LD). In addition, haplotypes containing HLA-DRB3~ HLA-DRB1*12:02 and HLA-C*07:06~ HLA-B*44:03~ HLA-DRB1*07:01 alleles were found as two independent protective factors for leprosy. Moreover, we replicated the previously associated HLA-DRB1*15:01 as leprosy risk factor and HLA-DRB1*04:05~HLA-DQA1*03:03 as protective alleles. When we narrowed the analysis to the single amino acid level, we found that the associations of the HLA alleles were largely captured by four independent amino acids at HLA-DRß1 positions 57 (D) and 13 (F), HLA-B position 63 (E) and HLA-A position 19 (K). Hence, analyses at the amino acid level circumvented the ambiguity caused by strong LD of leprosy susceptibility HLA alleles and identified four distinct leprosy susceptibility factors.

Family-based genome-wide association study of leprosy in Vietnam.

Leprosy is a chronic infectious disease of the skin and peripheral nerves with a strong genetic predisposition. Recent genome-wide approaches have identified numerous common variants associated with leprosy, almost all in the Chinese population. We conducted the first family-based genome-wide association study of leprosy in 622 affected offspring from Vietnam, followed by replication in an independent sample of 1181 leprosy cases and 668 controls of the same ethnic origin. The most significant results were observed within the HLA region, in which six SNPs displayed genome-wide significant associations, all of which were replicated in the independent case/control sample. We investigated the signal in the HLA region in more detail, by conducting a multivariate analysis on the case/control sample of 319 GWAS-suggestive HLA hits for which evidence for replication was obtained. We identified three independently associated SNPs, two located in the HLA class I region (rs1265048: OR = 0.69 [0.58-0.80], combined p-value = 5.53x10-11; and rs114598080: OR = 1.47 [1.46-1.48], combined p-value = 8.77x10-13), and one located in the HLA class II region (rs3187964 (OR = 1.67 [1.55-1.80], combined p-value = 8.35x10-16). We also validated two previously identified risk factors for leprosy: the missense variant rs3764147 in the LACC1 gene (OR = 1.52 [1.41-1.63], combined p-value = 5.06x10-14), and the intergenic variant rs6871626 located close to the IL12B gene (OR = 0.73 [0.61-0.84], combined p-value = 6.44x10-8). These results shed new light on the genetic control of leprosy, by dissecting the influence of HLA SNPs, and validating the independent role of two additional variants in a large Vietnamese sample.

Reducing Aspergillus fumigatus Virulence through Targeted Dysregulation of the Conidiation Pathway.

Inhalation of conidia of the opportunistic mold Aspergillus fumigatus by immunocompromised hosts can lead to invasive pulmonary disease. Inhaled conidia that escape immune defenses germinate to form filamentous hyphae that invade lung tissues. Conidiation rarely occurs during invasive infection of the human host, allowing the bulk of fungal energy to be directed toward vegetative growth. We hypothesized that forced induction of conidiation during infection can suppress A. fumigatus vegetative growth, impairing the ability of this organism to cause disease. To study the effects of conidiation pathway dysregulation on A. fumigatus virulence, a key transcriptional regulator of conidiation (brlA) was expressed under the control of a doxycycline-inducible promoter. Time- and dose-dependent brlA overexpression was observed in response to doxycycline both in vitro and in vivo. Exposure of the inducible brlA overexpression strain to low doses of doxycycline under vegetative growth conditions in vitro induced conidiation, whereas high doses arrested growth. Overexpression of brlA attenuated A. fumigatus virulence in both an invertebrate and mouse model of invasive aspergillosis. RNA sequencing studies and phenotypic analysis revealed that brlA overexpression results in altered cell signaling, amino acid, and carbohydrate metabolism, including a marked upregulation of trehalose biosynthesis and a downregulation in the biosynthesis of the polysaccharide virulence factor galactosaminogalactan. This proof of concept study demonstrates that activation of the conidiation pathway in A. fumigatus can reduce virulence and suggests that brlA-inducing small molecules may hold promise as a new class of therapeutics for A. fumigatus infection.IMPORTANCE The mold Aspergillus fumigatus reproduces by the production of airborne spores (conidia), a process termed conidiation. In immunocompromised individuals, inhaled A. fumigatus conidia can germinate and form filaments that penetrate and damage lung tissues; however, conidiation does not occur during invasive infection. In this study, we demonstrate that forced activation of conidiation in filaments of A. fumigatus can arrest their growth and impair the ability of this fungus to cause disease in both an insect and a mouse model of invasive infection. Activation of conidiation was linked to profound changes in A. fumigatus metabolism, including a shift away from the synthesis of polysaccharides required for cell wall structure and virulence in favor of carbohydrates used for energy storage and stress resistance. Collectively, these findings suggest that activation of the conidiation pathway may be a promising approach for the development of new agents to prevent or treat A. fumigatus infection.

Genetics of leprosy: today and beyond.

Leprosy is a chronic infectious disease of the skin and peripheral nerves that presents a strong link with the host genetic background. Different approaches in genetic studies have been applied to leprosy and today leprosy is among the infectious diseases with the greatest number of genetic risk variants identified. Several leprosy genes have been implicated in host immune response to pathogens and point to specific pathways that are relevant for host defense to infection. In addition, host genetic factors are also involved in the heterogeneity of leprosy clinical manifestations and in excessive inflammatory responses that occur in some leprosy patients. Finally, genetic studies in leprosy have provided strong evidence of pleiotropic effects between leprosy and other complex diseases, such as immune-mediated or neurodegenerative diseases. These findings not only impact on the field of leprosy and infectious diseases but also make leprosy a good model for the study of complex immune-mediated diseases. Here, we summarize recent genetic findings in leprosy susceptibility and discuss the overlap of the genetic control in leprosy with Parkinson's disease and inflammatory bowel disease. Moreover, some limitations, challenges, and potential new avenues for future genetics studies of leprosy are also discussed in this review.

Pleiotropic effects for Parkin and LRRK2 in leprosy type-1 reactions and Parkinson's disease.

Type-1 reactions (T1R) are pathological inflammatory episodes and main contributors to nerve damage in leprosy. Here, we evaluate the genewise enrichment of rare protein-altering variants in 7 genes where common variants were previously associated with T1R. We selected 474 Vietnamese leprosy patients of which 237 were T1R-affected and 237 were T1R-free matched controls. Genewise enrichment of nonsynonymous variants was tested with both kernel-based (sequence kernel association test [SKAT]) and burden methods. Of the 7 genes tested 2 showed statistical evidence of association with T1R. For the LRRK2 gene an enrichment of nonsynonymous variants was observed in T1R-free controls (PSKAT-O = 1.6 × 10-4). This genewise association was driven almost entirely by the gain-of-function variant R1628P (P = 0.004; odds ratio = 0.29). The second genewise association was found for the Parkin coding gene PRKN (formerly PARK2) where 7 rare variants were enriched in T1R-affected cases (PSKAT-O = 7.4 × 10-5). Mutations in both PRKN and LRRK2 are known causes of Parkinson's disease (PD). Hence, we evaluated to what extent such rare amino acid changes observed in T1R are shared with PD. We observed that amino acids in Parkin targeted by nonsynonymous T1R-risk mutations were also enriched for mutations implicated in PD (P = 1.5 × 10-4). Hence, neuroinflammation in PD and peripheral nerve damage due to inflammation in T1R share overlapping genetic control of pathogenicity.

Deciphering the genetic control of gene expression following Mycobacterium leprae antigen stimulation.

Leprosy is a human infectious disease caused by Mycobacterium leprae. A strong host genetic contribution to leprosy susceptibility is well established. However, the modulation of the transcriptional response to infection and the mechanism(s) of disease control are poorly understood. To address this gap in knowledge of leprosy pathogenicity, we conducted a genome-wide search for expression quantitative trait loci (eQTL) that are associated with transcript variation before and after stimulation with M. leprae sonicate in whole blood cells. We show that M. leprae antigen stimulation mainly triggered the upregulation of immune related genes and that a substantial proportion of the differential gene expression is genetically controlled. Indeed, using stringent criteria, we identified 318 genes displaying cis-eQTL at an FDR of 0.01, including 66 genes displaying response-eQTL (reQTL), i.e. cis-eQTL that showed significant evidence for interaction with the M. leprae stimulus. Such reQTL correspond to regulatory variations that affect the interaction between human whole blood cells and M. leprae sonicate and, thus, likely between the human host and M. leprae bacilli. We found that reQTL were significantly enriched among binding sites of transcription factors that are activated in response to infection, and that they were enriched among single nucleotide polymorphisms (SNPs) associated with susceptibility to leprosy per se and Type-I Reaction, and seven of them have been targeted by recent positive selection. Our study suggested that natural selection shaped our genomic diversity to face pathogen exposure including M. leprae infection.

Age-Dependent Association of TNFSF15/TNFSF8 Variants and Leprosy Type 1 Reaction.

A current major challenge in leprosy control is the prevention of permanent disabilities. Host pathological inflammatory responses termed type 1 reaction (T1R) are a leading cause of nerve damage for leprosy patients. The environmental or inherited factors that predispose leprosy cases to undergo T1R are not known. However, studies have shown an important contribution of host genetics for susceptibility to T1R. We have previously identified variants encompassing the TNFSF15/TNFSF8 genes as T1R risk factors in a Vietnamese sample and replicated this association in a Brazilian sample. However, we failed to validate in Brazilian patients the strong association of TNFSF15/TNFSF8 markers rs6478108 and rs7863183 with T1R that we had observed in Vietnamese patients. Here, we investigated if the lack of validation of these variants was due to age-dependent effects on association using four independent population samples, two from Brazil and two from Vietnam. In the combined analysis across the four samples, we observed a strong association of the TNFSF15/TNFSF8 variants rs6478108, rs7863183, and rs3181348 with T1R (pcombined = 1.5E-05, pcombined = 1.8E-05, and pcombined = 6.5E-06, respectively). However, the association of rs6478108 with T1R was more pronounced in leprosy cases under 30 years of age compared to the global sample [odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.54-2.46, pcombined = 2.5E-08 versus OR = 1.46, 95% CI = 1.23-1.73, pcombined = 1.5E-05]. A multivariable analysis indicated that the association of rs6478108 with T1R was independent of either rs7863183 or rs3181348. These three variants are known regulators of the TNFSF8 gene transcription level in multiple tissues. The age dependency of association of rs6478108 and T1R suggests that the genetic control of gene expression varies across the human life span.

A genome wide association study identifies a lncRna as risk factor for pathological inflammatory responses in leprosy.

Leprosy Type-1 Reactions (T1Rs) are pathological inflammatory responses that afflict a sub-group of leprosy patients and result in peripheral nerve damage. Here, we employed a family-based GWAS in 221 families with 229 T1R-affect offspring with stepwise replication to identify risk factors for T1R. We discovered, replicated and validated T1R-specific associations with SNPs located in chromosome region 10p21.2. Combined analysis across the three independent samples resulted in strong evidence of association of rs1875147 with T1R (p = 4.5x10-8; OR = 1.54, 95% CI = 1.32-1.80). The T1R-risk locus was restricted to a lncRNA-encoding genomic interval with rs1875147 being an eQTL for the lncRNA. Since a genetic overlap between leprosy and inflammatory bowel disease (IBD) has been detected, we evaluated if the shared genetic control could be traced to the T1R endophenotype. Employing the results of a recent IBD GWAS meta-analysis we found that 10.6% of IBD SNPs available in our dataset shared a common risk-allele with T1R (p = 2.4x10-4). This finding points to a substantial overlap in the genetic control of clinically diverse inflammatory disorders.

A genome wide association study identifies a IncRna as risk factor for pathological inflammatory responses in leprosy

Leprosy Type-1 Reactions (T1Rs) are pathological inflammatory responses that afflict a sub-group of leprosy patients and result in peripheral nerve damage. Here, we employed a family-based GWAS in 221 families with 229 T1R-affect offspring with stepwise replication to identify risk factors for T1R. We discovered, replicated and validated T1R-specific associations with SNPs located in chromosome region 10p21.2. Combined analysis across the three independent samples resulted in strong evidence of association of rs1875147 with T1R (p = 4.5x10-8; OR = 1.54, 95% CI = 1.32-1.80). The T1R-risk locus was restricted to a lncRNA-encoding genomic interval with rs1875147 being an eQTL for the lncRNA. Since a genetic overlap between leprosy and inflammatory bowel disease (IBD) has been detected, we evaluated if the shared genetic control could be traced to the T1R endophenotype. Employing the results of a recent IBD GWAS meta-analysis we found that 10.6% of IBD SNPs available in our dataset shared a common risk-allele with T1R (p = 2.4x10-4). This finding points to a substantial overlap in the genetic control of clinically diverse inflammatory disorders.

A Missense LRRK2 Variant Is a Risk Factor for Excessive Inflammatory Responses in Leprosy.

BACKGROUND: Depending on the epidemiological setting, a variable proportion of leprosy patients will suffer from excessive pro-inflammatory responses, termed type-1 reactions (T1R). The LRRK2 gene encodes a multi-functional protein that has been shown to modulate pro-inflammatory responses. Variants near the LRRK2 gene have been associated with leprosy in some but not in other studies. We hypothesized that LRRK2 was a T1R susceptibility gene and that inconsistent association results might reflect different proportions of patients with T1R in the different sample settings. Hence, we evaluated the association of LRRK2 variants with T1R susceptibility. METHODOLOGY: An association scan of the LRRK2 locus was performed using 156 single-nucleotide polymorphisms (SNPs). Evidence of association was evaluated in two family-based samples: A set of T1R-affected and a second set of T1R-free families. Only SNPs significant for T1R-affected families with significant evidence of heterogeneity relative to T1R-free families were considered T1R-specific. An expression quantitative trait locus (eQTL) analysis was applied to evaluate the impact of T1R-specific SNPs on LRRK2 gene transcriptional levels. PRINCIPAL FINDINGS: A total of 18 T1R-specific variants organized in four bins were detected. The core SNP capturing the T1R association was the LRRK2 missense variant M2397T (rs3761863) that affects LRRK2 protein turnover. Additionally, a bin of nine SNPs associated with T1R were eQTLs for LRRK2 in unstimulated whole blood cells but not after exposure to Mycobacterium leprae antigen. SIGNIFICANCE: The results support a preferential association of LRRK2 variants with T1R. LRRK2 involvement in T1R is likely due to a pathological pro-inflammatory loop modulated by LRRK2 availability. Interestingly, the M2397T variant was reported in association with Crohn's disease with the same risk allele as in T1R suggesting common inflammatory mechanism in these two distinct diseases.

Association of TNFSF8 regulatory variants with excessive inflammatory responses but not leprosy per se.

BACKGROUND: Type 1 reactions (T1R) affect a considerable proportion of patients with leprosy. In those with T1R, the host immune response pathologically overcompensates for the actual infectious threat, resulting in nerve damage and permanent disability. Based on the results of a genome-wide association study of leprosy per se, we investigated the TNFSF15 chromosomal region for a possible contribution to susceptibility to T1R. METHODS: We performed a high-resolution association scan of the TNFSF15 locus to evaluate the association with T1R in 2 geographically and ethnically distinct populations: a family-based sample from Vietnam and a case-control sample from Brazil, comprising a total of 1768 subjects. RESULTS: In the Vietnamese sample, 47 single-nucleotide polymorphisms (SNPs) overlapping TNFSF15 and the adjacent TNFSF8 gene were associated with T1R but not with leprosy. Of the 47 SNPs, 39 were cis-expression quantitative trait loci (cis-eQTL) for TNFSF8 including SNPs located within the TNFSF15 gene. In the Brazilian sample, 18 of these cis-eQTL SNPs overlapping the TNFSF8 gene were validated for association with T1R. CONCLUSIONS: Taken together, these results indicate TNFSF8 and not TNFSF15 as an important T1R susceptibility gene. Our data support the need for infection genetics to go beyond genes for pathogen control to explore genes involved in a commensurate host response.